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We apply the open systems concept and the influence 
functional formalism introduced in Paper I to establish a 
stochastic theory of relativistic moving spinless particles in a 
quantum scalar field. The stochastic regime resting between 
the quantum and semi-classical captures the statistical me- 
chanical attributes of the full theory. Applying the particle- 
centric world-line quantization formulation to the quantum 
field theory of scalar QED we derive a time-dependent (scalar) 
Abraham-Lorentz-Dirac (ALD) equation and show that it is 
the correct semiclassical limit for nonlinear particle-field sys- 
tems without the need of making the dipole or non-relativistic 
approximations. Progressing to the stochastic regime, we 
derive multiparticle time-dependent ALD-Langevin equations 
for nonlinear ly coupled particle-field systems. With these 
equations we show how to address time-dependent dissipa- 
tion/noise/renormalization in the semiclassical and stochastic 
limits of QED. We clarify the relation of radiation reaction, 
quantum dissipation and vacuum fluctuations and the role that 
initial conditions may play in producing non-Lorentz invariant 
noise. We emphasize the fundamental role of decoherence in 
reaching the semiclassical limit, which also suggests the correct 
way to think about the issues of runaway solutions and preac- 
celeration from the presence of third derivative terms in the 
ALD equation. We show that the semiclassical self-consistent 
solutions obtained in this way are "paradox" and pathology free 
both technically and conceptually. This self-consistent treat- 
ment serves as a new platform for investigations into problems 
related to relativistic moving charges. 



I. INTRODUCTION 

This is the second in a series of three papers |l]J|] ex- 
ploring the regime of stochastic behavior manifested by 
relativistic particles moving through quantum fields. This 
series highlights the intimate connections between back- 
reaction, dissipation, noise, and correlation, focusing on 
the necessity (and the consequences) of a self-consistent 
treatment of nonlincarly interacting quantum dynamical 
particle-field systems. 
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In Paper I we have set up the basic framework 
built on the influence functional and worldline quantiza- 
tion formalism. In this paper we apply the results obtained 
there to spinless relativistic moving particles in a quan- 
tum scalar field. The interaction is chosen to be the scalar 
analog of QED coupling so that we can avoid the com- 
plications of photon polarizations and gauge invariance in 
an electromagnetic field. The results here should describe 
correctly particle motion when spin and photon polariza- 
tion are unimportant, and when the particle is sufficiently 
decohered such that its quantum fluctuations effectively 
produce stochastic dynamics. 

The main results of this investigation are 

1 . First principle derivation of a time-dependent (modi- 
fied) Abraham-Lorentz-Dirac (ALD) equation as the 
consistent semiclassical limit for nonlinear particle- 
scalar-field systems without making the dipole or 
non-relativistic approximations. 

2. Consistent resolution of the paradoxes of the ALD 
equations, including the problems of runaway and 
acausal (e.g. pre-accelerating) solutions, and other 
"pathologies" . We show how the non-Markovian na- 
ture of the quantum particle open-system enforces 
causality in the equations of motion. We also discuss 
the crucial conceptual role that decoherence plays in 
understanding these problems. 

3. Derivation of multiparticle Abraham-Lorentz-Dirac- 
Langevin (ALDL) equations describing the quantum 
stochastic dynamics of relativistic particles. The 
familiar classical Abraham-Lorentz-Dirac (ALD) 
equation is reached as its noise-averaged form. The 
stochastic regime, characterized by balanced noise 
and dissipation, plays a crucial role in bridging the 
gap between quantum and (emergent) classical be- 
haviors. The N-particle-irreducible (NPI) and mas- 
ter effective action |?]-^| provide a route for gener- 
alizing our treatment to the self-consistent inclusion 
of higher order quantum corrections. 

We divide our introduction into three parts: First, a 
discussion of the pathologies and paradoxes of the ALD 
equation from the conventional approach and our sug- 
gested cures based on proper treatment of causal and non- 
Markovian behavior and self-consistent backreaction. Sec- 
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ond, misconceptions in the relation between radiation re- 
action, quantum dissipation and vacuum fluctuations and 
the misplaced nature and role of fluctuation-dissipation 
relations. Finally we give a brief summary of previous 
work and describe their shortcomings which justify a new 
approach as detailed here. 



There have been notable efforts to understand charged 
particle radiation-reaction in the classical and quantum 
theory (for a review see ||), of which extended charge 
theories ]l4j and quantum Langevin equations are impor- 
tant examples p5|^17[| (See Subsection C below). 



A. Pathologies, Paradoxes, Remedies and Resolutions 

The classical theory of moving charges interacting with 
a classical electromagnetic field has controversial difficul- 
ties associated with backreaction The generally ac- 
cepted classical equation of motion in a covariant form 
for charged, spinless point particles, including radiation- 
reaction, is the Abraham-Lorentz-Dirac (ALD) equation 

2" + (2e 2 /3m) {^z 2 + z") = {e/m)z v F^ t (z). (1.1) 

The timescale To = (2e 2 /3mc 3 ) determines the relative 
importance of the radiation-reaction term. For electrons, 
To ~ 10~ 24 sees, which is roughly the time it takes light 
to cross the electron classical radius ro ~ 10~ 15 m. The 
ALD equation has been derived in a variety of ways, often 
involving some regularization procedure that renormalizes 
the particle's mass 10 1. Feynman and Wheeler have de- 
rived this result from their "Absorber" theory which sym- 
metrically treats both advanced and retarded radiation on 
the same footing H . 

The ALD equation has strange features, whose status 
are still debated. Because it is a third order differential 
equation, it requires the specification of extra initial data 
(e.g. the initial acceleration) in addition to the usual posi- 
tion and velocity required by first order Hamiltonian sys- 
tems. This leads to the existence of runaway solutions. 
Physical (e.g. non -runaway) solutions may be enforced by 
transforming (1.1) to a second order integral equation with 
boundary condition such that the final energy of the par- 
ticle is finite and consistent with the total work done on it 
by external forces. But the removal of runaway solutions 
comes at a price, because solutions to the integral equa- 
tion exhibit the acausal phenomena of pre-acceleration on 
timescales r c . This is the source of lingering questions on 
whether the classical theory of point particles and fields is 
causalf]. 



B. Radiation Reaction and Vacuum Fluctuations 

The stochastic regime is characterized by competition 
between quantum and statistical processes: particularly, 
quantum correlations versus decoherence, and fluctuations 
versus dissipation. The close association of vacuum fluctu- 
ations with radiation-reaction is well-known, but the pre- 
cise relationship between them requires clarification. It is 
commonly asserted that vacuum fluctuations (VF) are bal- 
anced by radiation reaction (RR) through a fluctuation- 
dissipation relation (FDR). This can be misleading. Ra- 
diation reaction is a classical process (i.e., H = 0) while 
vacuum fluctuations are quantum in nature. A counter- 
example to the claim that there is a direct link between 
RR and VF is the uniformly accelerated charged parti- 
cle (UAP) for which the classical radiation reaction force 
vanishes, but vacuum fluctuations do not. However, at the 
stochastic level variations in the radiation reaction force 
away from the classical averaged value exist - it is this 
quantum dissipation effect which is related to vacuum fluc- 
tuations by a FDR. 

From discussions in Paper I we see how decoherence 
enters in the quantum to classical transition. Under rea- 
sonable physical conditions [18| the influence action or de- 
coherence functional is dominated by the classical solution 
(via the stationary phase approximation). Also from the 
' non-Mar kovian' fluctuation-dissipation relations found in 
Paper I, we see that the relationship between noise and dis- 
sipation is extremely complicated. Here too, decoherence 
plays a role in the emergence of the usual type fluctuation- 
dissipation kernel, just as it plays a role in the emergence 
of classicality. When there is sufficient decoherence, the 
FDR kernel (see Eq. (5.23) in Paper I) is dominated by 
the self-consistent semiclassical (decoherent) trajectories. 
An approximate, linear FDR relation may be obtained 
(see Eq. (5.24) in Paper I) except that the kernel is self- 
consistently determined in terms of the average (mean) 
system history, and the FDR describes the balance of fluc- 
tuations and dissipation about those mean trajectories. 

These observations have important ramifications on the 
appropriateness of assigning a FDR for radiation reaction 



1 Another mystery of the ALD equation is the problem of the 
constant force solution, where the charge uniformly accelerates 
without radiation reaction, despite it being well-established 
that there is radiation. This is a case where local notions of 
energy must be carefully considered in theories with local in- 
teractions (between point particles) mediated by fields ^2|Jl3) . 
The classical resolution of this problem involves recognizing 



that one must consider all parts of the total particle-field sys- 
tem with regard to energy conservation: particle energy, radia- 
tion energy, and non-radiant field energy such as resides in the 
so-called acceleration (or Shott) fields |13|. 
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and vacuum fluctuations. Ordinary (e.g. classical) radia- 
tion reaction is consistent with the mean particle trajecto- 
ries determined by including the backreaction force from 
the field. But at the classical level, there are no fluctua- 
tions, and hence no fluctuation-dissipation relation. Even 
regarding radiation reaction as a 'damping' mechanism is 
misleading- it is not 'dissipation' in the true statistical 
mechanical sense. Unlike a particle moving in some vis- 
cous medium whose velocity is damped, radiation reaction 
generally vanishes for inertial (e.g. constant velocity) par- 
ticles moving in vacuum fields (e.g. QED). Furthermore, 
radiation reaction can both 'damp ' and anti- 'damp ' parti- 
cle motion, though the average effect is usually a damping 
one. On the other hand, we show that the fluctuations 
around the mean-trajectory are damped, as described by 
a FDR. At the quantum stochastic level as discussed here, 
quantum dissipation effect results from the changes in the 
radiation reaction force that are associated with fluctu- 
ations in the particle trajectory around the mean, and 
therefore it is not the same force as the classical (i.e. av- 
erage) radiation reaction force. Only these quantum pro- 
cesses obey a FDR. This is a subtle noteworthy point. 

C. Prior work in relation to ours 

f. There are many works on nonrelativistic quantum 
(and semiclassical) radiation reaction (for atoms as 
well as charges moving in quantum fields) including 
those by Rohrlich, Moniz-Sharp, Cohen- Tannoudji 
et. al., Milonni, and others fl|l3| , |l9| -pT| . Kampen, 
Moniz, Sharp, and others mM suggested that the 
problems of causality and runaways can be resolved 
in both classical and quantum theory by considering 
extended charge models. However, this approach, 
while quite interesting, misses the point that point 
particles (like local quantum field theories) obey a 
good low energy effective theory in their own right. 
While it is true that extended objects can cure in- 
finities (e.g., extended charge models, string the- 
ory), it is nonetheless important to recognize that 
at sufficiently low energy QED as an effective field 
theory is 'consistent' without recourse to its high- 
energy limit. Wilson, Weinberg and others have 
shown how effective theory description is sufficient 
to understand low energy physics because compli- 
cated, and often irrelevant, high energy details of the 
fine structure are not being probed at the physical 
energy scales of interest. The effective theory ap- 
proach has also provided a new perspective on the 
physical meaning of renormalization and the source 
of divergences, showing why nonrenormalizable in- 
teractions are not a disaster, why the apparent large 
shift in bare-parameters resulting from divergent 
loop diagrams (i.e. renormalization) doesn't inval- 
idate perturbation theory, and clarifying how and 



when high-energy structure does effect low-energy 
physics. Therefore one should not need to invoke 
extended charge theory to understand low energy 
particle dynamics. 

Our work is a relativistic treatment which includes 
causal Non-Markovian behavior and self-consistent 
backreaction. We think this is a better approach 
even in a nonrelativistic context because the regular- 
ized relativistic theory never 'breaks down'. In our 
approach we adopt the ideas and methods of effec- 
tive field theory with proper treatment of backreac- 
tion to account for the effects of high-energy (short- 
distance) structure on low-energy behavior, and to 
demonstrate the self-consistency of the semiclassical 
particle dynamics. We consider the consequences 
of coarse-graining the irrelevant (environmental) de- 
grees of freedom, and the nature of vacuum fluc- 
tuations and quantum dissipation in the radiation 
reaction problem. We emphasize the crucial role of 
decoherence due to noise in resolving the pathologies 
and paradoxes, and in seeing how causal QED leads 
to a (short-time modified) causal ALD semiclassical 
limit. 

2. By examining the time-dependence of operator 
canonical commutation relations Milonni showed the 
necessity of electromagnetic field vacuum fluctua- 
tions for radiating nonrelativistic charges 

13- The 

conservation of the canonical commutation relations 
is a fundamental requirement of the quantum theory. 
Yet, if a quantum particle is coupled to a classical 
electromagnetic field, radiative losses (dissipation) 
lead to a contraction in "phase-space" for the ex- 
pected values of the particle position and momen- 
tum, which violates the commutation relations, ft is 
the vacuum field that balances the dissipation effect, 
preserving the commutation relations [x, p] = ih 
as a consequence of a fluctuation-dissipation rela- 
tion (FDR). The particle-centric worldline frame- 
work allows us to consider similar issues for rela- 
tivistic particles. Interestingly, our covariant frame- 
work shows how quantum fluctuations (manifesting 
as noise) appear in both the time and space coor- 
dinates of a particle. This is not surprising since 
relativity requires that a physical particle satisfies 
= t 2 (r) — z 2 (t) = f , and therefore any vari- 
ation in spatial velocities must be balanced by a 
change in i° keeping the particle "on-shell" . 

3. The derivation and use of quantum Langevin equa- 
tions (QLE's) to describe fluctuations of a system 
in contact with a quantum environment has a long 
history. Typically, QLE's are assumed to describe 
fluctuations in the linear response regime for a sys- 
tem around equilibrium, but its validity does not 
need be so restrictive. Nonequilibrium conditions 
can be treated with the Feynman- Vernon influence 
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functional. Caldeira and Leggett's study of quantum 
Brownian motion (QBM) has led to an extensive lit- 
erature |23|] , particularly in regard to decoherence 
issues [|24[| . Barone and Caldeira [p5] have applied 
this method to the question of whether nonrelativis- 
tic, dipole coupled electrons decohere in a quantum 
electromagnetic field. An advantage of Barone and 
Caldeira's work is that it is not limited to initially 
factorized states; they use the preparation function 
method which allows the inclusion of initial par- 
ticle field correlations. Despite this, Romero and 
Paz have pointed out that the preparation function 
method still suffers from an (implicit) unphysical de- 
piction of instantaneous measurement characterizing 
the initial state preparation |26| . 

4. Ford, Lewis, and O'Connell have extensively dis- 
cussed the electromagnetic field as a thermal bath 
in the linear, dipole coupled regime |l7|, and pio- 
neered the application of QLE's to nonrelativistic 
particle motion in QED. They have detailed the con- 
ditions for causality in the thermodynamic, equilib- 
rium limit described by the late-time linear quantum 
Langevin equation. A crucial point of their analy- 
sis is that particle motion can be (depending on the 
cutoff of the field spectral den sity ) runaway free and 
causal in the late- time limit ^7J. In |2qj , they sug- 
gest a form of the equations of motion that gives 
fluctuations without dissipation for a free electron, 
but this result is special to the particular choice of a 
field cutoff made by the authors, namely, one imply- 
ing that the bare mass of the particle exactly van- 
ishes. When one does not assume a special value 
for the cutoff one generally does find both fluc- 
tuations and dissipation, though Ford and Lewis's 
counter-example shows how careful one must be in 
handling FDRs. Another counter-example against 
making hasty generalizations is the vanishing of the 
averaged (classical) radiation reaction for uniformly 
accelerated particles despite quantum fluctuations 
that make the (suitably coarse-grained) trajectories 
stochastic. 

We do agree with the result that the particle motion 
is only causal and consistent when the field cutoff is 
below a certain critical valued] (determined from the 
particles classical radius) , but we find no compelling 
justification for the claim that the cutoff should take 
exactly that special value. 



2 A noted example is in Caldeira and Leggett's original work 
Jlt| . Their master equation governs a reduced density matrix 
which is not positive definite because they choose an (infinite) 
cutoff greater than the inverse temperature (see Hu, Paz, and 
Zhang [g3| for details). 



It is certainly interesting that with the special field- 
cutoff chosen by Ford and O'Connell higher time- 
derivatives vanish from the equations of motion. 
However, even without such a specially chosen cut- 
off, the influence of higher derivative terms are 
strongly suppressed at low energies. On this issue 
we take the effective theory point of view which em- 
phasizes the generic insensitivity of low-energy phe- 
nomena to unobserved high-energy structures. Since 
higher-time derivatives in the equations of motion 
do not (necessarily) violate causality, there is no 
reason to rule them out. It is perhaps instructive 
to review the meaning of renormalizable field theo- 
ries in the light of effective theories. Effective the- 
ories generically have higher derivative interactions 
that are again strongly suppressed at low energies. 
What makes renormalizable terms special is that 
these are the interactions that remain relevant (or 
marginal) at low energies, while non-renormalizable 
terms (such as higher derivative terms) are expo- 
nentially suppressed, hence, it is no longer believed 
that non-renormalizable terms are fundamentally 
excluded so long as one recognizes that physical 
theories are mostly effective theories. Since QED 
is certainly an effective theory, these observations 
(together with the suppression of the higher time- 
derivative terms) weaken the argument that the spe- 
cial cutoff theory of FOL is fundamental. Ford and 
O'Connell also propose a relativistic generalization 
of their modified equations of motion for the aver- 
age trajectory derived from the nonrelativistic QLE 
[ p9[ . Our derivation of the stochastic limit from a 
relativistic quantum mechanics (i.e. worldline for- 
mulation) and field theory goes well beyond this by 
allowing the treatment of fully quantum relativis- 
tic processes, and by deriving a relativistic Langevin 
equation. 

5. A perturbative expansion of the ALD equation up 
to order e 3 has been derived from QED field the- 
ory by Krivitskii and Tsytovich p0[ , including the 
additional forces arising from particle spin. Their 
work shows that the ALD equation may also be un- 
derstood from field theory, but the authors have not 
addressed the role of fluctuations, correlation, deco- 
herence, time-dependent renormalization, nor self- 
consistent backreaction. Our derivation yields the 
full ALD equation (all orders in e for the perturba- 
tion expansion employed by J|(|) through the loop 
(semiclassical) expansion; this method automati- 
cally includes all tree-level diagrammatic effects. 

6. Low pi] ] showed that runaway solutions apparently 
do not occur in spin 1/2 QED. But Low does not 
derive the ALD equation, nor address the semiclassi- 
cal/stochastic limit. We also emphasize, again, that 
it does not, and should not, matter whether parti- 
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cles are spin 1/2, spin 0, or have some other internal 
structure in regard to the causality of the low energy 
effective theory for center of mass particle motion. 

7. Using the influence functional, Diosi |32| derives 
a Markovian master equation in non-relativistic 
quantum mechanics. In contrast, it is our intent 
to emphasize the non-Markovian and noncquilib- 
rium regimes with special attention paid to self- 
consistency. The work of |f32|| differs from ours in the 
treatment of the influence functional as a functional 
of particle trajectories in the relativistic worldline 
quantization framework. Ford has considered the 
loss of electron coherence from vacuum fluctuation 
induced noise with the same noise kernel that we 
employ |33| ]. However, his application concerns the 
case of fixed or predetermined trajectories. 

D. Organization 

In Section II we obtain the influence functional (IF) for 
spinless relativistic particles. We then define the stochas- 
tic effective action (S s t oc h) for this model, using our re- 
sults from Paper I, and use it to derive nonlinear integral 
Langevin equations for multiparticle spacetime motion. In 
Section III we consider the single particle case. We derive 
the (scalar-field) Abraham-Lorentz-Dirac (ALD) equation 
as the self-consistent semiclassical limit, and show how this 
limit emerges free of all pathologies. In Section IV we de- 
rive a Langevin equation for the stochastic fluctuations of 
the particle spacetime coordinates about the semiclassical 
limit for both the one particle and multiparticle cases. We 
also show how a stochastic version of the Ward-Takahashi 
identities preserves the mass-shell constraint. In Section 
V, we give a simple example of these equations for a single 
free particle in a scalar field. We find that in this particu- 
lar case the quantum field induced noise vanishes, though 
this result is special to the overly-simple scalar field, and 
does not hold for the electromagnetic field ^ . In the final 
section, we summarize our main results and mention areas 
of applications, of both theoretical and practical interest. 

II. SPINLESS RELATIVISTIC PARTICLES 
MOVING IN A SCALAR FIELD 

A. The model 

Relativistic quantum theories are usually focused pri- 
marily on quantized fields, the notion of particles following 
trajectories is somewhat secondary. As explained in Paper 
I B, we employ a "hybrid" model in which the environ- 
ment is a field, but the system is the collection of particle 



spacetime coordinates z% (r„) (i.e., worldlines) where n in- 
dicates the n th particle coordinate, with worldline param- 
eter t„, charge e n , and mass m n . The free particle action 
is 

S M = Yl I dr nV / SfK + V(z n ( Tn ))} . (2.1) 

From Sa follow the relativistic equations of motion: 

mz£/y/£* = -d»V(z n ). (2.2) 

For generality, we include a possible background potential 
V, in addition to the quantum scalar field environment. 
The scalar current is 

j[z,x) = X! / dT n e n u n (T n )S(x - z„(r„)), (2.3) 
where 



u n (r n ) = yi« (T„)i„, M (r n ). (2.4) 

In Paper III [^|, we treat a vector current coupled to the 
electromagnetic vector potential A^. In both cases we 
assume spinless particles. The inclusion of spin or color 
is important to making full use of these methods in QED 
and QCD. 

In the full quantum theory reparametrization invari- 
ance of the relativistic particle plays a crucial role, and 
the particle-field interaction must respect this symmetry. 
When dealing with the path integral for the quantized 
particle worldline, a quadratic form of the action is often 
more convenient than the square- root action in ( |2.l| ). The 
worldline quantum theory is a gauge theory because of the 
constraint that follows from reparametrization invariance 
under r — > t' (t) = r + e (r) , making this perhaps the 
simplest physical example of a generally covariant theory. 
(We may therefore view this work as a toy model for gen- 
eral relativity, and for string theory, in the semiclassical 
and stochastic regime.) The path integral must therefore 
be gauge-fixed to prevent summing over gauge-equivalent 
histories. In |34| we treat the full quantum theory, in de- 
tail, and derive the "worldline influence functional" after 
fixing the gauge to the so-called proper-time gauge. This 
give the worldline path integral in terms of a sum over 
particle trajectories (r) satisfying z 2 = N 2 , but where 
the path integral still involves an integration over possible 
N (the range of this integration depends on the bound- 
ary conditions for the problem at hand). However, apply- 
ing the semiclassical (i.e., loop) expansion in j34j] we find 
that the N — 1 trajectory gives the stationary phase solu- 
tion; hence, the semiclassical limit^ corresponds to setting 



3 Note however that in the stochastic limit there can poten- 
tially be fluctuations in the value of the constraint. 
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N = 1. In the following we will therefore assume the con- 
straint z 2 — -/V 2 , or even z — 1, when we are discussing 
the semiclassical limit for trajectories. We show that this 
is self-consistent in that it is preserved by the equations 
of motion including noise-induced fluctuations in the par- 
ticle trajectory. Since we are addressing the semiclassi- 
cal/stochastic limit in this paper there is no need to adopt 
the quadratic form of action for a relativistic particle; it 
is adequate, and it turns out to be simpler in this p aper , 
to just work with the scalar functions u n defined in ( |2.4| ), 
and the square-root form of the relativistic particle action 
in Q. 

While we emphasize the microscopic quantum origins, 
we might also view our model in analogy to the treat- 
ment of quantum fields in curved spacetime. There, one 
takes the gravitational field (spacetime) as a classical sys- 
tem coupled to quantum fields. One important class of 
problems is then the backreaction of quantum fields on 
the classical spacetime. The backreaction of their mean 
yields the semiclassical Einstein equation which forms the 
basis of semiclassical gravity. The inclusion of fluctuations 
of the stress-energy of the quantum fields and the induced 
metric fluctuations yields the Einstein-Langevin equations 
|55) which forms the basis of this stochastic semiclassical 
gravity |36|. In our work here, the particle coordinates 



are analogous to the gravitational field (metric tensor). 
Whereas the Einstein-Langevin equations have not been 
derived from first principles for the lack of a quantum the- 
ory of gravity, we shall take advantage of the existence of 
a full quantum theory of particles and fields to examine 
how, and when, a regime of stochastic behavior emerges. 



B. The influence functional and stochastic effective 
action 

The interaction between the particles and a scalar field 
is given by the monopole coupling term 

Sint = J dxj[z,x)ip{x) (2.5) 

= e 5Z / dT n ll n (T n )ip(z n (T n )). 

This is the general type of interaction treated in Paper 
I. In the second line we have used the expression for the 
current (^3). The expression for the influence functional 
(Eq. (3.16) of Paper I) is, 

=exp{ -^[dxf dx'[2j-(x)G R (x,x')j + (x') 
-ij-(x)G H (x',x')j-(x')]}, (2.6) 



where 



3 = -j[z',x)) 
j+ = {j[z,x)+j[z',x))/2 



(2.7) 



and {z, z'} are a pair of particle histories. G R and G H 
are the scalar field retarded and Hadamard Green's func- 



tions, respectively. Substitution of (2.3) then gives the 
multiparticle influence functional 

F[{z} , {z'}\ = exp{ - e — [ dr n f dr m (2.8) 

x6(T-z° n )8(z°-z° m ) 

X [^i/ n G (^n; z m)^m ^rS~* (^n> z m)^m 
— U n (j (z n , z rn )u m -+- U n (j \z n -, z m )u m 
~t~ Un,G (z n , Z m )i/ m (^n> z m)^m 

+ u n G (z n , z rn )u m ~ u n G (z n , z to )m„J ^ 

where u' n = u n (z' n ). The influence functional may be ex- 
pressed more compactly by using a matrix notation where 



(u n ,u' n ) = (u^ul), giving 
F[z a ] = cxp{ dT n dT m 

nm 

x K^nm U m + u n^nm u m) } 

= 6XP I ~H ( U " G ™™ Um + U n G nmUr, 



exp<j-S 7F [z a ] 



(2.9) 

(2.10) 
(2.11) 



The superscript T denotes the transpose of the column 
vector u. In ( p. 10 ), and below, we leave the sum S„ m and 
integrations Jrfr^dr m implicit for brevity. The matrices 
G" m ,G R m are given by 



G 



H 



e(T-z n )e(z n -z° m ) 

(Gfn) { z n' Z rn) — ^Yl2) i Z n' Z rn) 
~ ^(21) i Z n' Z rn) ^(22) ( Z n > Z rn) 



(2.12) 



and 



G R m = 9(T-z n )e(z° n -z° m ) 

1^(11) i Z m z rn) ^(12) { Z n> z m) 
~ ^(21) \ Z n' Z m) ~ ^(22) ( Z n) z m) 



(2.13) 



In Eq. (2.11), we have defined the influence action Sif. 
qH,r are -j-^g scalar-field Hadamard/Retarded Green's 
functions evaluated at various combinations of spacetime 



points 



The stochastic effective action is defined by (see Jy) 
S x [z±] = SaIz*} + J dxj-^^Xxix) (2.14) 
+ 2 f dx'G R {x,x')j + [z ± ,x')) 



U n G Jlm U m 



(i 



where 



U^Wn = U n (z n ) X (Z n ) ~ U n (z' n ) X (4) 



(2.15) 



and x ( z n) is a (classical) stochastic field (see Section 3 of 
Paper I) evaluated at the spacetime position of the nth 
particle. The stochastic field has vanishing mean and au- 
tocorrelation function given by 

(x (*) X (x')) 8 = K({<P (x) , (x')}) = hG H (x, x') . (2.16) 

Hence, x i x Y s statistics encodes those of the quantum 
field tp (x) . 

The stochastic effective action provides an efficient 
means for deriving self-consistent Langevin equations that 
describe the effects of quantum fluctuations as stochastic 
particle motion for sufficiently decohered trajectories (for 
justification in using S s t oc h in this way see Paper I). 



C. Langevin integrodifferential equations of motion 

Deriving a specific Langevin equation from the general 
formalism of Paper I requires evaluating the cumulants 
Cq , defined in Eq. (4.11) of [Q. Because the multiparti- 
cle case involves a straightforward generalization, we begin 
with the single particle theory. The noise cumulants are 
then defined as 



g 



s q s IF [f 



V 5z& (r n ) ...Szm (r m ) 



(2.17) 



=o 



with Sif given in (2.11). The superscript (?y) indicates 
that these are cumulants for an expansion of Sjp with 
respect to the stochastic variables rj^ (t) . 

We define new variables v = Uu where v T =(u 1 ,w 2 ) 
and u T = (u, u') , and 



u=if 1 1 
2 V 1 1 



Then the influence action has the form 

^[z Q ] = ^(v T Gfv + zv T Gfv) 



(2.18) 



(2.19) 



where 
G R = 



UG H U T (2.20) 

/ r<R n R j.n R n R n R ^r< R j-n R 

I <-ri 2 -t- ( - r 21 — U 22 ( - r ll+ ( - T 12 + ( - T 21 + Lr 22 

\ r< R n R n R j^c r n R j^c r n R n R 

\ Wl - Lt 12 — ( - t 21+ Lt 22 ( - z ll+ ( - r 12 — ( - r 21 — ( - r 22 



and 

G? = UG ff U T 



(2.21) 



GH 1 s~*H 1 s~iH 1 /~iH /~tH r^H i 

11 ' l - T 12+ LT 21+ l - r 22 U ll~ ( - T 12~r t - Z 21 — U 22 

GH 1 /~iH /~iH /^iH /~*H r~iH /~<H , /~iH 

ll+ l - T 12 — ( - T 21 — l - T 22 "-^ll - Lr 12 — Lr 21+ ( - r 22 



The lowest order cumulant in the Langevin equation, 
C[^, is found by evaluating SSif/Sz^~ , and then set- 
ting z~ — 0. There are two kinds of terms that arise: 
those where 5/8z^~ acts on v, and those where it acts on 
G R ' . For linear theories, G R ' H is not a function of the 
dynamical variables, but is instead (at most) a function 
of some predetermined kinematical variables that a pri- 
ori specify the trajectory. Hence, there is no contribution 
from SG R ' H / 'Sz^~ . For the <5 v/(5z M ~ te rms, setting z~ = 
collapses the matrices ( 2.20| and 2.21 ) to one term each: 
only the (1,2) term of G R , and the (1,1) term of G R , 
survive. But, setting z~ — gives v 1 = 0, v = u, and 
since the (1,1) term of G R is proportional to two factors 
of v 1 , it also vanishes. When 5/5z^~ acts on G R , only 
the (2,2) terms survive (because it is the only term not 
proportional to a factor of v 1 ). For similar reasons, the 
only contributing element of G H is also its (2, 2) term. To 

evaluate 



(«Z? (M) /& 



=0) we note 



6 

Tz 



5 

Jz~^' 



5z->* 

After a little algebra, we are left with 

SG r (z(t),z(t')) 



/Sz^ 



(2.22) 



(2.23) 



where the derivative only acts on the z(t), and not the 
z(t') argument, in G R . The same algebra in evaluating 

the ( $Gv(22)/^ z>l ~^ \z-=o term shows that all the factors 

cancel, and therefore G H does not contribute to the first 
cumulant at all. Because the imaginary part of Sif[z°"] 
does not contribute to the first cumulant, the equations of 
motion of the mean-trajectory are explicitly real, which is 
an important consequence of using an initial value formu- 
lation like the influence functional (or closed-time-path) 
method. 

The first cumulant, describing radiation reaction, is 
therefore given by 



SS 



IF 



8z^-(t) 



(2.24) 



f T dr'e 2 { 


( Su(z) \ 




\6z^ (r) J 



G r {z(t),z{t'))u{z(t')) 



■u{z{t)) 



5G r (z(t),z(t')) 
5zV( t ) 



i(z{t'))}. 



This expression is explicitly causal both because the 
proper-time integration is only over values r' < r, and 
because of the explicit occurrence of the retarded Green's 
function. Contrary to common perception, the radiation 
reaction force given by C$ is not necessarily dissipative 

in nature. For instance, we shall see that vanishes 
for uniformly accelerated motion, despite the presence of 
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radiation from a uniformly accelerating charge. In other 
circumstances, [z, t) may actually provide an anti- 
damping force for some portions of the particle trajectory. 

Hence, radiation reaction is not a purely dissipative, 
or damping, process. The relationship of radiation reac- 
tion to vacuum fluctuations is also not simply given by 
a fluctuation-dissipation relation^. In actuality, only that 
part of the first cumulant describing deviations from the 
purely classical radiation reaction force is balanced by fluc- 
tuations in the quantum field via generalized fluctuation 
dissipation relations. This part we term quantum dissi- 
pation because it is of quantum origin and different from 
classical radiation reaction. Making this distinction clear 
is important to dispel common misconceptions about ra- 
diation reaction being always balanced by vacuum fluctu- 
ations. 

Next we evaluate the second cumulant. After similar 
manipulations as above, we find that the second cumulant 
involves only G H . We note that the action of Sj/Sz^ on 
an arbitrary function is given by 



8z^(t) 
dx I dr 



(2.25) 



Sz^ 



{[z 2 ) 1/2 8{x-z{r)))f{x) 



J±(J !L _ 



u(t) 



Because of the constraint z 2 — N 2 , it follows that u = 0, 
and we may set u — N. Also, = 0. Then 



N 



{'z^ + i' y i^<9„]} 



f{z{r)) 

(2.26) 



This last expression defines the operator (z) . We have 
used d/dr = z v d v . The operator satisfies the identity 



ifvtn (z) = z" 



N 



— {z li + z"z [ll d l/] } 



(z> i z ll + z»z u z [fJ ,d v] ) =0. 



(2.27) 



This identity ensures that neither radiation reaction nor 
noise- induced fluctuations in the particle's trajectory move 



4 Fluctuation-dissipation related arguments have been fre- 
quently invoked for radiation reaction problems, but these in- 
stances generally entail linearizing (say, making the dipole ap- 
proximation), and/or the presence of a binding potential that 
provide a restoring force to the charge motion making it peri- 
odic (or quasi-periodic). Radiation reaction in atomic systems 
is one example of this the assumption of some kind of 

harmonic binding potential provides another example M,Bgfl. 



the particle off-shell (i.e., the stochastic equations of mo- 
tion preserve the constraint z 2 — N 2 for any constant TV). 
With these definitions, the noise ^(t) is defined by 



77 m (t) = eh 1/2 w^{z) X (z(T)) 



zh 1 ' 2 



(2.28) 



and the second-order noise correlator by 



C^[z;r,T')^({^(r),^(r')}) (2.29) 
= e 2 hw^z)w»( Z i)({ X (z(T)) ,x(z(r'))}) 
= e 2 hw»(z)w l, {z')G H {z(T),z{T 1 )). 

The operator (z) acts only on the z in G H (z, z') ; like- 
wise, the operator w v (z f ) acts only on z' . 

This scalar field result is reminiscent of electromag- 
netism, where the Lorentz force from the (antisymmetric) 
field strength tensor F^ v is = z^F^ = z u d^A v ]. The 
antisymmetry of F^ v implies z^T^ = 0. We may define 
a scalar analog of the antisymmetric (second rank) field 
strength tensor by 



px 

fj,h> 



~~[» 



d v ]X- 



(2.30) 



This shows that the second term on the right hand side of 
( 2.28 ) gives the scalar analog of (a stochastic) electromag- 



netic Lorentz force: 



z v F*. The first term on the 



RHS of ( 2.28 ), z^x ( z ) > does not occur in the treatment of 
the electromagnetic field. In the scalar-field theory, this 
term may be thought of as a stochastic component to the 
particle mass. 

The stochastic equations of motion are 



mz^ = -d^V{z) + e 2 J dT'w^z)G R (z(T), z(t')) 

+ eh 1 / 2 w ll {z) X {z{T)). (2.31) 



The result (2.31) is formally a set of nonlinear stochastic 
integrodifferential equations for the particle trajectories 
Zu[t); t). Noise is absent in the classical limit found by the 
prescription H — > Q (this definition of classicality is for- 
mal in that the true semiclassical/classical limit requires 
coarse-graining and decoherence, and is not just a matter 
of taking the limit % — * 0). 



The generalization of ( 2.31 ) to multiparticles is now 
straightforward. If we re-insert the particle number in- 
dices, the first cumulant is 



Cl x [z,t) 



SS IF [z] 



z m =0 



6z£-(t) 
= / &r' m <?w nll {z n ) 

X G^(z„(t„), Z m (T m ))u(z m (T m ), 



(2.32) 



S 



the noise term is given by 

<(r) = eh^w^xizj^)), (2.33) 
and the noise correlator is 

C^: m) [z;r n ,rU) (2.34) 

= (K(^),c(r;)}) 

= e 2 hw^{z n )w U m {z m )G H {z n (T n ), Z m {T' m )). 

The nonlinear multiparticle Langevin equations are there- 
fore 

mzn^ (t) (2.35) 
= -d^V(z n (r)) + eS 1 / 2 w M (z n ) X (2„(T)) 



dT'wJz n {T))G R (z n (r),z m (T')). 



= 2 E 



The n ^ m terms in ( p.3q ) are particle-particle inter- 
action terms. Because of the appearance of the retarded 
Green's function, all of these interactions are causal. The 
n = m terms are the self-interaction (radiation-reaction) 
forces. The n ^ m noise correlator terms represent nonlo- 
cal particle-particle correlations: the noise that one parti- 
cle sees is correlated with the noise that every other parti- 
cle sees. The nonlocally correlated stochastic field x( x ) re- 
flects the correlated nature of the quantum vacuum. From 
the fluctuation-dissipation relations found in Paper I (Eq. 
(5.23) and Eq. (5.24)), the n = m quantum noise is re- 
lated to the n = m quantum dissipative forces. Under 
some, but not all, circumstances^] the n ^ m correlation 
terms are likewise related to the n ^ m propagation (in- 
teraction) terms through a multiparticle generalization of 
the FDR, called a propagation-correlation relation [P9| . 



We have already noted that the first term in (2.28) has 
the form of a stochastic contribution to the particle's effec- 
tive mass, thus allowing us to define the stochastic mass 



as 



m v z^ = 



m + eh}-l 2 x{z) ) z 



(2.36) 



Fluctuations of the stochastic mass automatically preserve 
the mass-shell condition since (m-^z^z^) \ z=i — 0. Like- 



wise, the effective stochastic force F* satisfies 



0. (2.37) 



5 See and Paper I, Section IV, for a discussion of this point. 
In brief, the noise correlation between spacelike separated 
charges does not vanish owing to the nonlocality of quantum 
theory, but the causal force terms involving G R does always 
vanish between spacelike seperated points. These two kinds 
of terms are only connected through a propagation/correlation 
relation (the multiparticle generalization of an FDR) when one 
particle is in the other's casual future (or past). 



These results show that the stochastic fluctuation-forces 
preserve the constraint i 2 = N 2 . 



III. THE SEMICLASSICAL REGIME: THE 
SCALAR-FIELD ALD EQUATION 



We have emphasi zed i n Pape r 1 that the emergence of a 
Langevin equation (2.31 or 2.35) presupposes decoherence, 



which works to suppress large fluctuations away from the 
mean-trajectories. In the semiclassical limit, the noise- 
average vanishes, and the semiclassical equations of mo- 
tion for the single particle are given by the mean of (2.31). 
The stochastic regime admits noise induced by quantum 
fluctuations around these decohered mean solutions. In 
our second series of papers we explore the stochastic be- 
havior due to higher-order quantum effects that refurbish 
the particle's quantum nature. 

On a conceptual level, we note that the full quan- 
tum theory in the path integral formulation involves sum- 
ming over all worldlines of the particle^] joining the ini- 
tial and final spacetime positions, Zi and z/, respectively. 
Furthermore, the generating functional for the worldline- 
coordinate expectation values (see Subsection 3. A of Paper 
I jjj) involve a sum over final particle positions Zf. In these 
path integrals, there is no distinction between, say, run- 
away trajectories and any other type of trajectoryn. Fur- 
thermore, no meaningful sense of causality is associated 
with individual (i.e. fine-grained or skeletonized) histories 
in the path integral sum. Any particular path in the sum 
going through the intermediate point z (r) at worldlinc 
parameter time r bears no causal relation to it than going 
through the point z'(r') at some later parameter time r'. 

Where questions of causality, uniqueness, and runaways 
do arise is in regard to the solutions to the equations of mo- 
tion for correlation function of the worldline-coordinates, 
particularly for the expectation value giving the mean tra- 
jectory. It is here that an initial value formulation for 
quantum physics is crucial because only then are the equa- 
tions of motion guaranteed to be real and causal. In con- 
trast, equations of motion found from the in-out effective 
action (a transition amplitude formulation) arc generally 
neither real nor casual ^|0]]. Moreover, the equations of 



6 Actually, just what kinds of histories/ worldlines are allowed 
in the sum is determined by both the gauge choice in the world- 
line quantization method, and by the boundary conditions. 

7 In fact, trajectories in the path integral can be even stranger 
than the runaways that appear in classical theory. Most paths 
are non-differentiable (infinitely rough) and may even include 
those that travel outside the lightcone and backward-in-time 
(see previous footnote). Clearly, non-uniqueness of paths is 
not an issue in the path integral context. 
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motion for correlation functions must be unique and fully 
determined by the initial state if the theory is complete. 

These observations lead to a reframing of the questions 
that are appropriate in addressing the semiclassical limit 
of quantum particle-field interactions. Namely, what are 
the salient features of the decoherent, coarse-grained histo- 
ries which qualify as semiclassical particle motion? Quan- 
tum theory permits the occurrence of events which would 
be considered classically improbable or forbidden in the 
particle motions (e.g., runaways, tunneling, or apparently 
acausal behavior). These effects are allowed because of 
quantum fluctuations. Bearing in mind the requirement 
of classicality, we need to show that out of the infinite pos- 
sibilities in the quantum domain the interaction between 
particle and field (including the self-field of the particle) 
is causal in the observable semiclassical limit. 

So the pertinent questions in a quantum to stochastic 
treatment are the following: What are the equations of 
motion for the mean and higher-order correlation func- 
tions? Are these equations of motion causal and well- 
defined? How significant are the quantum fluctuations 
around the quantum-averaged trajectory? When does de- 
coherence suppress the probability of observing large fluc- 
tuations in the motion? When do the quantum fluctua- 
tions assume a classical stochastic behavior? It is only by 
addressing these questions that the true semiclassical mo- 
tion may be identified, together with the noise associated 
with quantum fluctuations which is instrumental for deco- 
herence. With this discussion as our guide, we proceed in 
two steps. First, in this section, we find the semiclassical 
limit for the equations of motion. Second, in the follow- 
ing section, we describe the stochastic fluctuations around 
that limit. Since we arc dealing with a nonlinear theory, 
the fluctuations themselves must depend on the semiclas- 
sical limit in a self-consistent fashion. These two steps 
constitute the full backreaction problem for nonlinear par- 
ticle field interactions in the semiclassical and stochastic 
regimes. For more general conceptual discussions on deco- 
herent histories and semiclassical domains, see |l^ , ^ , [l3| . 



A. Divergences, regularized QED, and the initial state 

Even in classical field theory, Green's functions are usu- 
ally divergent, which calls for regularization. Ultraviolet 
divergences arise because the field contains modes of ar- 
bitrarily short wavelengths to which a point particle can 
couple. These facts have contributed to considerable con- 
fusion about the correct form the radiation reaction force 
can take on for a charged particle. There are also in- 
frared divergences that arise from the artificial distinction 
between soft and virtual photon emissions. These diver- 
gences are an artifact of neglecting recoil on the particle 
motion. 

In this work, we take the effective field theory philos- 
ophy as our guide. One does not need to know the de- 



tailed structure of the correct high-energy theory because 
low-energy processes are largely insensitive to them 22 
beyond the effective renormalization of system parame- 
ters and suppressed high-energy corrections. Take QED 
as an example, where relevant physics occurs at the physi- 
cal energy scale E p h ys ■ In QED calculations, it is generally 
necessary to assume a high-energy cutoff E cuto f f that reg- 
ularizes the theory. One also assumes a high-energy scale 
E e ff well above E p h ys where presumably new physics is 
possible, thus E p h ys « E e ff < E cuto ff- The effective 
theory that describes physics at scales well below E e ff is 
found by coarse-graining (integrating out) the high-energy 
fluctuations between E e ff and E cuto ff- These fluctuations 
renormalize the masses and charges in the effective theory, 
and add new local, effective-interactions that reproduce 
the effects of the coarse-grained physics between E e f f and 
Ecutoff - These so-called nonrenormalizablc effective inter- 
actions are suppressed by some power of (E p h ys / E e f f), 
making them usually insignificant compared to the renor- 
malizable terms in the effective theory. This justifies the 
viability of ordinary QED, which is comprised of just 
those renormalizable terms. It further ensures that physics 
based on the physical masses and charges probed at E p h ys 
will be insensitive to modifications to QED above E e ff. 

Following the standard procedure in QED, we now pro- 
ceed to regularize the Green's functions by modifying its 
high-frequency components. Two traditional methods of 
regularization are the dimensional and Pauli-Villars tech- 
niques, but neither of these are fully adequate for the prob- 
lem at hand. Dimensional regularization is usually applied 
after an analytic continuation to the Euclidean path inte- 
gral, but here we work with the real-time path integral and 
nonequilibrium boundary conditions. It turns out that the 
Pauli-Villars method does not suppress high-energy field 
modes strongly enough for our purposes because it poten- 
tially allows high-energy corrections to modify the long- 
time particle motion. 

Ultimately, there is no reason to prefer any one regu- 
larization scheme over another beyond issues of suitabil- 
ity and particularities (or for preserving some particular 
features or symmetries of the high-energy theory). For 
our purposes, a sufficiently strongly regulated retarded 
Green's function is given by the Gaussian smeared func- 
tion 



G%(a)=O(z o (T)-z°(T'))0(a) 



(3.1) 



with support inside and on the future lightcone of z (r') 
We have defined 



a (s) = (s) y p (s) 
y» (s) ee z» (r) - z» (r>) 



(3.2) 
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In the limit A — ► oo, 



lim G% (a) = 5 (a) /2tt, 

A-^oo 



(3.3) 



giving back the unregulated Green's function. Expanding 
the function a for small s gives 



cr (s) = l/ (s) Vll (s) 



(3.4) 



and 




(3.5) 



assuming a timelike trajectory z^ (r) . 

Now that we have defined a cutoff version of QED, we 
briefly reconsider the type of initial state assumed in this 
series of papers. The exact influence functional is generally 
found for initially uncorrelated particle(s) and field, where 
the initial density matrix is 



P (U) = Pz (U) ® Ph> 



(3.6) 



Such a product state, constructed out of the basis |z,t) ® 
\tp(x)), is highly excited with respect to the interacting 
theory's true ground state. If the theory were not ultra- 
violet regulated, producing such a state would require in- 
finite energy. This is not surprising since it is impossible 
to take a fully dressed particle state and strip away all 
arbitrarily-high-energy co rrela tions by some finite energy 
operation. With a cutoff, ( |3.6| ) is now a finite energy state 
and thus physically achievable, though depending on the 
cutoff A it may still require considerable energy to pre- 
pare. The important point, which we show below, is that 
the late-time behavior is largely insensitive to the details 
of the initial state, while the early-time behavior is fully 
causal and unique. 



B. Single-particle scalar-ALD equation 

The semiclassical limit is given by the equations of mo- 
tion for the mean-trajecto ry. This is just the solution 



to the Langevin equation (2.31) without the noise term 
where, as we indicated earlier, t is the classical proper- 
time such that z 2 = A^ 2 . The regulated semiclassical equa- 
tions of motion are then 



m z M (r) - d,V(z) = C£ } (3.7) 
= e 2 f dr' (z, (r) + j? (r) z { , (r) d v] ) Gf{z{r), z(r')) 

where the <9„ acts only on the z (t) in Gj±, and not the 
z (r') term. We add the subscript to too because it is 



the particle's bare mass. This nonlocal integrodifferen- 
tial equation is manifestly causal by virtue of the re- 
tarded Green's function. Eq. (3.7) describes the full 



non-Markovian semiclassical particle dynamics, which de- 
pend on the past particle history for proper-times in the 
range [r, n] . These second order integrodifferential equa- 
tions have unique solutions determined by the ordinary 
Newtonian (i.e. position and velocity) initial data at ti. 

T herefore, the nonlocal equations of motion found from 
(3/7) do not have the problems associated with the classi- 
cal Abrah am-L orentz-Dirac equation. But, as an integral 
equation, ( |3.7| ) depends on the entire history of z (t) be- 
tween Ti and r. For this reason, the transformation from 



(3.7) to a local differential equation of motion requires 
every derivative of z (t) (e.g., d n z(r) /dr n ). This is the 
origin of higher-time-derivatives in the local equations of 
motion for radiation-reaction (classically as well as semi- 
classically) . By integrating out the field variables, we have 
removed an infinite set of nonlocal (field) degrees of free- 
dom in favor of nonlocal kernels whose Taylor expansions 
give higher-derivative terms. 



We may now transform the integral equation (3.7) to 
a local differential equation by expanding the functions 
y^ (s) around s = 0. In Figure 1, we show a semiclassical 
trajectory with respect to the lightcone at z (r) . Taking r 
as fixed, we change integration variables using ds = dr' . 
Next, we need the expansions 



-yi> (s) = (t) s + F (r) s 2 /2 
+ z» (r) s 3 /6 + O (s 4 ) , 

-if (s) = z" (t) + z" (t) s 

+ z» (r),s 2 /2 + G( S 3 ), 

^f^(s) = s - z 2 s 3 /24 + o (s 4 ) 



(3.8) 



(3.9) 



We set z 2 = 1 for convenience since we are concerned 
with the semiclassical limit. Also, we use = 0, and 

z 2 = — z^ z n to simplify as necessary. We denote both 
dz/dr by z and dy/ds by y. Using 

da (s) 



dz^ 
da (s) 
ds 



2»"(a) (3.10) 
2y"y ll = 2 S -z- 2 s 3 /2 + 0(s 4 ) 



allows us to write the gradient operator as 
da d 



<>,, =- - = 2y 



dz^ da M V ds J ds 



da\ 1 d 



(3.11) 



y v y u ds ' 



We also define 



r = t - n 



(3.12) 



to be the total elapsed proper-time for the particle since 
the initial time n. Recall that Ti is defined by z° (n) = ti 
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where ti is the initial time at which the state p(U) = 
Pz (ti) ® Pip (U) is defined. 

These definitions and relations let us write the RHS of 
0)as 



da(uPG2(8) + uP~G%(8) (3.13) 



M 2ds 
-G«( S ) + ....). 



6 c?s 



where 



(n + 1)! ds 



< 3 =^ (r) 

< 2) = (^ 2 + *m) 



«W=4 n) (i,z,...,d n+1 z/dr" +1 ) 

The functions depend on higher-derivatives of the par- 
ticle coordinates up through order n + 1. We next define 
the time-dependent parameters 



h(r) 



dsGl{s) = ah(l- 



r(l/4,A 4 r 4 /2) 

iWi) 



(3.15) 



and 



ff (n) (r) 



rfs "7 77 -V 

(n + l)!ds 



(3.16) 



32 -(n-a)/4 - n , 
— - _ 7 (l + -LV/2), 



7r 3 / 2 (n + l) 

where T (x) is the Gamma function, T (x, y) is the incom- 
plete Gamma function, and 7 (x, y) = T (x) — T(x,y) . The 
constant a = 31" (5/4) / (2 1 / 4 tt 3 / 2 ) ~ .41 depends on the 
details of the high-energy cutoff, but is of order one for 
reasonable cutoffs. 

The coefficients are bounded for all r, including the limit 
r — > 00. In fact, 



lim g {n) (r) = lim g (n) (r) = 0, 



(3.17) 



taking either limit. Therefore, we are free to exchange the 
sum over n and the integration over s to find a convergent 
series expansion of (3.7). The resulting local equations of 
motion are 



m z„ (r) - d^V(z) = C {r,) 



(3.18) 



e 2 Mr) U W+ e 2 5>(™) (r) U W. 



From (3.14) we see that the terms give time-dependent 
mass renormalization; we accordingly define the renormal- 
ized mass as 

m (r) = mo — e 2 h (r) — e 2 g^ (r) 

too + 5m (r) . (3.19) 

Similarly, the term is the usual third derivative ra- 
diation reaction force from the ALD equation. Thus, the 
equations of motion may be written as 



m (r) Z/j, (r) - d M V(z) 



(r) 



(3.20) 



where 



(3.14) (r) 



(r) (z, t z 2 + Zp) + e 2 J2 



l (n) (r)«W 



(3.21) 



These semiclassical equations of motion are almost of 
the Abraham-Lorentz-Dirac form, except for the time- 
dependent renormalization of the effective mass, the time- 
dependent (non-Markovian) radiation reaction force, and 
the presence of higher than third-time-derivative terms. 
The mass renormalization is to be expected, indeed it oc- 
curs in the classical derivation as well. What an initial- 
value formulation reveals is that renormalization effects 
are time-dependent owing to the nonequilibrium dynam- 
ics of particle-field interactions. The reason why renor- 
malization is viewed as a time-independent procedure is 
because one usually only deals with equilibrium conditions 
or asymptotic in-out scattering problems. 



IV. RENORMALIZATION, CAUSALITY, AND 
RUNAWAYS 

A. Time-dependent mass renormalization 

The time-dependence of the effective mass to (r) (and 
particularly the radiation reaction coefficients g' n ') (r)) 
play crucial roles in the semiclassical behavior. The ini- 
tial particle at time ti (r = 0) is assumed, by our choice 
of initial state, to be fully uncorrelated with the field. 
This means the particle's initial mass is just its bare mass 
mo. Interactions between the particle and field then re- 
dress the particle state, one of the consequences being 
that the particle acquires a cloud of virtual field excita- 
tions (quanta) that change its effective mass. For a par- 
ticle coupled to a scalar field there are two types of mass 
renormalization effects, one coming from the z^tp interac- 
tion, and the other coming from the z^z^d^ip interaction. 
We commented in Section II. C that the latter interaction 
is essentially a scalar field version of the electromagnetic 
coupling, if one defines the scalar field analog of the field 
strength tensor 



= Z[udn]<P- 
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We find that this interaction contributes a mass shift at 
late times of 



„a„ ]V = ae 2 A/2. 



(4.1) 



In Paper III it is shown that Sm^a^y is just one half 
the mass shift that is found when the particle moves in 
the electromagnetic field. This is consistent with the in- 
terpretation of the electromagnetic field as equivalent to 
two scalar fields (one for each polarization); therefore one 
expects, and finds, twice the mass renormalization in the 
electromagnetic case. (For the same reason we see below 
that the radiation reaction force, and stochastic noise, are 
each also reduced by half compared to the electromagnetic 
case.) 

But the scalar field, unlike the electromagnetic field, 
also has a z^ip interaction that gives a mass shift Srriz V 



-25 



The combined mass shift is then negative. 



At late-times, the renormalized mass is 
TOqo = lim to (r) 

r — >oc 

= to + SrriEM + Sm sca i ar 
= too + ae 2 A/2 - aAe 2 
= toq — ae 2 A/2. 



(4.2) 



For a fixed mo, if the cutoff satisfies A > 2mo/ae 2 , the 
renormalized mass becomes negative and the equations of 
motion become unstable. It is well-known that environ- 
ments with overly large cutoffs can qualitative change the 
system dynamics. In the following we will assume that the 
bare mass and cutoff are such that TOqq > 0, which in turn 
implies that too > and gives a bound on the cutoff: 

A < 2m /ae 2 

(which is of the order of the inverse classical electron ra- 
dius). We emphasize that for the scalar field case the field 
interaction reduces the effective particle mass, whereas in 
the electromagnetic (vector) field case the field interaction 
increases the effective particle mass (see paper III 0]). If 
we followed the FOL p7| , |27| -p9t suggestion of picking a cut- 
off A = 2too /ae 2 we would then conclude that TOqo = as 
a consequence of the opposite sign in the mass renormal- 
ization effect. This is a kind of "critical" case balanced 
between stable and unstable particle motion, and while 
there may be special instances where such behavior is of 
interest, it does not seem to represent the generic behavior 
of typical (massive) particles interacting with scalar field 
degrees of freedom. We will defer further comparison with 
EM quantum Langevin equations to Paper III Q where 
the mass renormalization has the more familiar sign. 

We should further elaborate, however, on the role of 
causality in the equations of motion. In the preceding 
paragraph we see that for the long-term motion to be sta- 
ble (i.e. positive effective particle mass) there is an upper 
bound on the field cutoff. This is well-known, as detailed 



in |17| , |27| -|29|, for example. There, the QLE is assumed to 
be of the form 

tux -f f dt'fx (t - t') x (t') + V' (x) =F(t), (4.3) 



where the time integration runs from — oo to t. Causality 
is related to the analyticity of jl (z) (the Fourier trans- 
form of (t — t')) in the upper half plane o f the complex 
variable z. The traditional QLE of the type ( |4.3| ) with the 
lower integration limit set as U = — oo effectively makes 
the assumption that the particle has been in contact with 
its environment far longer than the bath memory time. In 
our analysis the lower limit in (3.7) is Tj, not — oo, because 



it is the (nonequilibrium) dynamics at early times immedi- 
ately after the initial particle/field state is prepared where 
we have something new to say about how causality is pre- 
served, and runaways are avoided. 

The proper-time dependence of the mass renormaliza- 
tion is shown in Figure 2. The horizontal axis marks the 
proper time in units of the cutoff timescale 1 /A. The verti- 
cal scale is arbitrary, depending on the particle mass. No- 
tice one unusual feature: the mass shift is not monotonic 
with time, but instead first overshoots its final asymp- 
totic value. This occurs because there are two competing 
mass renormalizing interactions that have slightly differ- 
ent timescales. The time-dependent mass-shift is a rapid 
effect, the final dressed mass is reached within a few 1/A. 



B. Nonequilibrium radiation reaction 

The constants g n (r) determine how quickly the particle 
is able to rebuild its own self-field, which in turn controls 
the backreaction on the particle motion, ensuring that it is 
causal. We have found, as is typical of effective field the- 
ories, that the equations of motion involve higher deriva- 
tive terms (i.e. d n z/dr n for n > 3) beyond the usual ALD 
form. In fact, all higher derivative terms are in principle 
present so long as they respect the fundamental symme- 
tries (e.g., Lorentz and reparametrization invariance) of 
our particle-field model. The timescales and relative con- 
tributions of these higher-derivative forces are determined 
by the coefficients g n (r) . For a fixed cutoff A, the late- 
time behavior of the g n (r) scale as 

2"/ 2 r (i + n/2) A2 _ n 

9" ( r ) = '3/2 A ( 4 - 4 ) 

7 — ► oo (271") 

Therefore, the g<i term has the late-time limit 

92(r) = l, (4.5) 

r — ► oo ^ n 

which is independent of A. In the terminology of effective 
field theory, this makes the n = 2 term renormalizable 
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and implies that it corresponds to a 'marginal' coupling 
or interaction. In the same terminology the mass renor- 
malization, scaling as A, corresponds to a relevant (renor- 
malizable) coupling. In this vein, the higher-derivative 
corrections proportional to g n >2, suppressed by powers of 
the high-energy scale A, correspond to so-called irrelevant 
couplings. 

If we had written the most general action for 
the particle-field system consistent with Lorentz and 
reparametrization invariance, we should have included 
from the beginning higher-derivative interactions. The ef- 
fective field theory methodology when carried out fully 
shows that these high-derivative interactions are in fact 
irrelevant at low-energies. Therefore, the ordinary ALD 
form for radiation reaction (e.g., the n — 2 term) results 
not as a consequence of the details of the 'true' high-energy 
theory, but because any Lorentz and reparametrization 
invariant effective theory behaves like a renormalizable 
particle-field theory that has the ALD equation as the 
low-energy limit. With this we come to the important 
conclusion that for a wide class of theories with the ap- 
propriate symmetries the low-energy limit of radiation re- 
action for (center of mass) particle motion is given by the 
ALD result. For theories with additional structure (such 
as spin), there will of course be additional radiation reac- 
tion terms important at the energy scale associated with 
the additional structure. 

In Figure 3 we show the proper-time dependence of 
(/ 2 ) (r), with the same horizontal timescale as in Fig- 
ure 2. It is evident that the radiation reaction force ap- 
proaches the ALD form quickly, essentially on the cutoff 
time. While the radiation reaction force is non-Markovian, 
the memory time of the field-environment is very short. 
The non-Markovian evolution is characterized by the de- 
pendence of the coefficients m (r) and (r) on the ini- 
tial time Ti , but become effectively independent of r, after 
t~a ~ 1/A. This is an example of the well-known behav- 
ior for quantum Brownian motion models found in p4| . 
The effectively transient nature of the short-time behavior 
helps justify our use of an initially uncorrelated (factoriz- 
able) particle and field state. 

Also notice the role played by the particle's elapsed 
proper-time, r = t — Ti. For particles with small aver- 
age velocities^ the particle's proper-time is roughly the 
background Minkowski time coordinate t. In this case, 
the time-dependences of the renormalization effects are 
approximately the same as one would find in the non- 
relativistic situation, since r = (r — Tj) ~ (t — ti). 
For rapidly moving particles time-dilation can signifi- 



The small velocity approximation is relative to a choice of 
reference frame. In paper I (see Appenix A) we discuss how 
the choice of initially factorized state at some time ti picks out 
a special frame. 



cantly lengthen the observed timescales with respect to 
the background Minkowski time. This indicates that 
highly relativistic particles will take longer (with respect 
to background time t) to equilibrate with the quantum 
field/environment than do more slowly moving particles. 
Finally, we note that the radiation reaction force is exactly 
half that found for the electromagnetic field. 



C. Causality and early-time behavior 

We have shown that at low-energies and late-times the 
ALD radiation reaction force is generic (neglecting other 
particle structure, such as spin, which gives rise to ad- 
ditional low-energy corrections). But the classical ALD 
equation is plagued with pathologies like acausal and run- 
away solutions. The equations of motion in the form ( [3.7| ) 
are clearly unique and causal, they require only ordinary 
initial data, and do not have runaway solutions. However, 
it is i nstructive to see how the equations of motions in the 
form (3.20 3.21), involving higher derivatives, preserve the 
causal nature of the solutions with radiation reaction. To 
address these questions we now examine the early-time 
behavior of the equations of motion. 

The presence of higher derivatives in ( |3.21 ) is an in- 
evitable consequence of converting from a set of nonlocal 
integr al equation (3.7) to a set of local differe ntial equa- 
tions ( |3.20| , |3.2lD . As differential equations, ( |3.20| , |3.2l[ ) 
may seem to be unphysical at first-sight because of 



the apparent need to specify initial data 



» 



(Ti) 



d n z fl (ri) /d,T n for n > 2. But the coupling constants g n (r) 
satisfy the crucial property that 



g n in) = 0, 



(4.6) 



for all n. We see this graphicall y for <?2 (?*) in Figure 3, 
and analytically for all g n (r) in ( 3.16| ). Consequently, the 
particle self-force at t = identically vanishes (to all 
orders), and only smoothly rebuilds as the particle's self- 
field is reconstituted. Therefore, the initial data at Ti is 
fully (and uniquely) determined by 



mz^ (n) = -9 M V (z (n)) , 



(4.7) 



which only requires the ordinary Newtonian initial data. 
The initial values for the higher derivative terms (e.g. n > 
2) in (3.21) are determined iteratively from 



1 dr? 



d n - 



dT?' 



d»V(z(n)) 



-2 M 



(4.8) 



Given the (Newtonian) initial data, the equations of mo- 
tion are determined uniquely for all later times, to any 
order in n. In the classical ALD equations, one finds run- 
aways even in the case of vanishing external potential, 
V = 0. In our case, because g n (n) = for all n, if V = 0, 
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Z/F (tj) = for all n > 2. With these initial conditions 



(and V = 0), the equations of motion (3.20, 3^2]] ) arc 



with unique solutions 



Zfj, (r) = z M (t;) + (t - Ti) (n) . 



(4.9) 



(4-10) 



Runaways do not arise. 

As in the classical derivation of the Abraham-Lorentz- 
Dirac equation, these semiclassical equations of motion are 
equivalent to a second order differential-integral equation. 
In the classical case the integral equations are "derived" 
from the third-order ALD equation. If one chooses bound- 
ary conditions for the integral equation that eliminate run- 
aways, one inevitably introduces pre-accelerated solutions. 
In our c ase, we d erive the higher-derivative equations of 
motion ( 3.2C , |3.21 ) from the explicitly casual non-local in- 
tegral equation (3.7). In doing this we avoid introducing 



unphysical solutions at short times as a consequence of 
the initially vanishing time-dependent parameters g n (r) . 
This turns out to be true for any good regularization of 
the field's Green functions. 



V. THE STOCHASTIC REGIME AND THE 
ALD-LANGEVIN EQUATION 

A. Single particle stochastic limit 



The nonlinear Langevin equation in (2.31) shows a com- 
plex relationship between noise and radiation reaction. A 
nonlinear Langevin equation similar to these have been 
found in stochastic gravity by Hu and Matacz |35|] de- 
scribing the stochastic behavior of the gravitational field 
in response to quantum fluctuations of the stress-energy 
tensor. We now apply the results in Paper I (see Sec- 
tion IV §) giving the linearized Langevin equation for 
fluctuations around the mean trajectory. In the Langevin 



equation (2.31) we define z M = z M + z 1 *, where z 11 is the 



semiclassical solution from the previous section, and we 
work to first order in the fluctuation variable zf". The free 
kinetic term is given by 

dr'r (/) ( P Sz [2] 



\5z^{t)5z v {t') 



m(r) / dT l z v {T')g flv —5{T-T') 



m (r) 



d 2 Zfj, (r) 
dr 2 ' 



(5.1) 



where we have included the time-dependent mass renor- 
malization effect in the kinetic term. The external poten- 
tial term is 

[«rrt(J^)-W*™, (5.2, 



6z»(t)5z v (t') 



dz» 2 



hence, the second derivative of V acts as a force linearly 
coupled to z. The dissipative term for z involves the (func - 
tional) derivative of the first cumulant (see Eq. ( |3.18 )) 
with the mass renormalization (n = 1) term removed, as it 
has already been included in the kinetic term ( |5.1[ ) above. 
We therefore need the (functional) derivative with respect 



to Zfj, of the radiation reaction term 



R.R. 



r 



To O (5) , we have 



E 



g (n) ( r ) 



,(«) 



. , „ (*) 

dTZAT) ^zVW 



(5.3) 



(5.4) 



We note that f R R (z) is a function of the derivatives z^ 
with n > 1, therefore we have 



oo a fR , R .(z\ 



n=l 
m+l 



dz y : 



{n) 



(5.5) 



e 2 E 4 n) W £ 



<?( m ) (r) 94 m) 



m=2 



dzi n) 



where df R R /dz^ are partial derivatives with respect to 
the n th derivative of z (r) . The sum over n terminates at 
m+l since only contains derivatives up through order 
z (™+ 1 )_ Hence, the Langevin equations for the fluctuations 
z (around z) are 



m R (t) 



d% (r) 



dr 2 "» 

m+l oo 

»7m W + e 2 E E 

n—1 m— 2 



S 2 V(z) 

Sz^ 2 
'gW (r) du { ™ y 
ml dzl n) 



(5.6) 



(r) 



These are solved with z found from ( 3.20 ). The noise 
?7 M (r) is given by 



Vn ( t ) = w^xiz (r)) 



(z^ (r)+ z"z [v d^x(z(T)), 

(5.7) 



where the lowest order noise is evaluated using the semi- 
classical solution z. 



The correlator (r) r\ v (r')) (see Eq. 5.9) is then found 
using the field correlator along the mean-trajectory: 

({ X (z(t)),x (z (/))}) = hG H (z (r) , z (/)) . (5.8) 



Equation (5.7) shows that the noise experienced by the 
particle depends not just on the stochastic properties of 
the quantum field, but on the mean-solution z (r) . For 



instance, the z^ x( z ) term in (5.7) gives noise that is pro- 
portional to the average particle acceleration. The second 
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term in (5.7) depends on the antisymmetrized combina- 
tion fr^ d v ]x(z). It follows immediately from (x) = that 
(»7(t))=0. 

From (|5.9[), we see that the noise correlator is 



(»7p ( T ) V» ( T ')> 

= ^(T)C(r'){x(^(r))x(2(T'))) 

= (e 2 h/c 4 ) ^ + I % <9f]J (X + 9*j 

x G ff (z,z') 



:,, + | A ^iHi: + r'^4] (5.9) 



' ~ V 2 da 
x (47rA c r c )G ff (a). 



y 2 da 



We have defined f/^ = z M (r) = z M (r') and a = y 2 . The 
noise (when the field is initially in its vacuum state) is 
highly nonlocal, reflecting the highly correlated nature of 
the quantum vacuum. Only in the high temperature limit 
does the noise become approximately local. Notice that 
the noise correlator inherits an implicit dependence on the 
initial conditions at Tj through the time-dependent equa- 
tions of motion for z. But when tA ^> 1, the equations of 
motion become effe ctiv ely independent of the initial time. 
The noise given by (5.9) is not stationary except for special 
cases of the solution z(r). 

Since the n > 2 terms represent high-energy corrections 
suppressed by powers of A, the low-energy Langevin equa- 
tions are well approximated by keeping only the n = 2 
term. The linear dissipation (n = 2) term becomes, using 



,(2) 



._ „2 

(z M z + Zn), 



e 2 £ 4 B) W 



n=l 

_2 



V 2 > (r) du$ Y 
4 2! d£ 



(n) 



(5.10) 



+ z 



where 



g {2} (r) (z z 2 fgp„- ZpZ^j 

(gnu- zp,zSj) 

9 V (r) ( o a h P0 f T) , 



(5.11) 



and a 2 —z . This lowest order expression for dissipation 
in the equations of motion for z^ involves time derivatives 
of z through third order. However, just as was the case 
for the semiclassical equations of motion, the vanishing 
of g( 2 ^ (r) for r = implies that the initial data for the 
equations of motion are just the ordinary kind involving 
no higher than first order derivatives. The lowest-order 



late-time (Ar >> 1) single particle Langevin equations 
arc thus 



rj (t) = m z^ (t) + z v 



_ dV(z) 



dz^dz v 



(5.12) 



B. Multiparticle stochastic limit and stochastic Ward 
identities 

It is a straightforward generalization to construct multi- 
particle Langevin equations. The two additional features 
are particle-particle interactions, and particle-particle cor- 
relations. The semiclassical limit is modified by the addi- 
tion of the terms 



dT m w^ [z n (r„))Gf (z„(t„) ))] . (5.13) 



The use of the regulated Ga is essential for consistency 
between the radiation that is emitted during the regime 
of time-dependent renormalization and radiation-reaction; 
it ensures agreement between the work done by radiation- 
reaction and the radiant energy. A field cutoff A implies 
that the radiation wave front emitted at Tj is smoothed on 
a time scale A -1 . 

We note one significant difference between the single 
particle and multiparticle theories. For a single parti- 
cle, the dissipation is local (when the field is massless), 
and the semiclassical limit (obtaining when rA ^> 1) is 
essentially Markovian. The multiparticle theory is non- 
Markovian even in the semiclassical limit because of mul- 
tiparticle interactions. Particle A may indirectly depend 
on its own past state of motion through the emission of 
radiation which interacts with another particle B, which in 
turn, emits radiation that re-influences particle A at some 
point in the future. The nonlocal field degrees of freedom 
stores information (in the form of radiation) about the 
particle's past. Only for a single timelike particle in flat 
space without boundaries is this information permanently 
lost insofar as the particle's future motion is concerned. 
These well-known facts make the multiparticle behavior 
extremely complicated. 



We may evaluate the integral in (5.13) using the identity 



6 ({z n - z m ) 2 ) = 5(t - t*)/2 (z n - z m ) u z mVj (5.14) 

where r* is the time when particle m crosses particle n's 
past lightcone. (The regulated Green's function would 
spread this out over approximately a time A -1 ). Defin- 
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mg z nm = z n - z m , we use 



dT rn d m G (z n . ) 

dG R (z nm ) 



(5.15) 



dTr? 



1 d 

2 (-^nm) ^mf 



+ / dsG — 

as 



(-2-nm) ^tj 



The particle-particle interaction terms are then given by 



[^ + z z^d^G R (z n ,z m ) 



(5.16) 



2 (^nm) 2 m( y C?T n 



where, as before, t he (z n ) acts only on the z n , and not 
the z m (t*) . Eq. ( 5.16| ) are just the scalar analog of the 
Licnard-Wiechert forces. They include both the near-field 
and far-field effects. 

The long range particle -part icle terms in the Langevin 
equations are found using ( 2.35 ) , so we have the additional 
Langevin term for the nth particle, found from 



e 2 J2 J 

h i2) A 

cy f— \ Cx — Aq- 



(z ) X 



(5.17) 



\Zrvmj Zrnft 



- ' 'Tm—T* 



Equation (5.17) contains no third (or higher) derivative 
terms. The multiparticle noise correlator is 

(Un ( T n) Vm (jm)}) = < ( T «) ™m ( T ™) ({X (*n) X (^m)}) • 

(5.18) 

In general, particle-particle correlations between spacelike 
separated points will not vanish as a consequence of the 
nonlocal correlations implicitly encoded by G . For in- 
stance, when there are two oppositely charged particles 
which never enter each other's causal future, there will 
be no particle-particle interactions mediated by G R ; how- 
ever, the particles will still be correlated through the field 
vacuum via G . This shows that it will not be possible 
to find a generalized multiparticle fluctuation-dissipation 
relation under all circumstances J39|. 

We now briefly consider general properties of the 
stochastic equations of motion. Notice that the noise sat- 
isfies the identity 

z n ( T ) Vnn (t) = % (t) W nU {z n )x{z n ) = 0, (5.19) 



which follows as a consequence of z£w nfl = 0, for each 
particle number separately. This is an essential property 
since the part icle fluctuations are real, not virtual. We 
may use fl5.19| ) to prove what might be called (by analogy 
with QED) stochastic Ward- Takahashi identities Q. The 
n-point correlation functions for the particle-noise are 



({Vnx(n)-V^i (n) -Vu„ ( T n)}) 

= UJ W (Zi) ...W Ut (Zi) ...W Un (z n ) 
X ({X(zi) -X(Zz) -X(Zn)}) , 



(5.20) 



with each w Ui (zi) acting on the corresponding x( z i)- 
These correlation functions may both involve different 
times along the worldline of one particle (i.e. self-particle 
noise) , and correlations between different particles. From 
(KM and (pi), follow "Ward" identities: 



for all i and n. 

(5.21) 



The contraction of on-shell momenta (i M ) with the 
stochastic correlation functions always vanishes. These 
identities are fundamental to the consistency of the rela- 
tivistic Langevin equations. 



VI. EXAMPLE: FREE PARTICLES IN THE 
SCALAR FIELD VACUUM 

As a concrete example, we find the Langevin equations 
for V{z) = and the scalar field initially in the vacuum 
state. Using the stochastic equations of motion we can 
address the question of whether a free particle will expe- 
rience Brownian motion induced by the vacuum fluctua- 
tions of the scalar field. When V(z) — 0, we immedi- 
ately see that for the semiclassical equations, z (r) = 
constant; the radiation- reaction force identically vanishes 
(for all orders of n). The fact that g n (0) = uniquely 
fixes the constant as zero due to the initial da ta b ound- 
ary conditions for the higher derivative terms (|4.8j); and 
hence, there are no runaway solutions characterized by 
a non-zero constant. We may write z (r) = and 
z^ (r) = , where is a constant spacetime velocity 
vector which satisfies the (mass-shell) constraint v 2 = 1. 
A particle moving in accordance with the semiclassical so- 
lution neither radiates nor experiences radiation-reactionR. 



Using a 



(r — r') , we immediately find 



9 In constast, the mean-solution for a free particle in a quan- 
tum scalar field found in [^(J does not possess the Galilean in- 
variance of our results, and its motion is furthermore damped 
proportionally to its velocity. 
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from ( 2.31 ) and (5.6) (using a 2 —z = 0) the (linear-order) 
Langevin equation 



m 



+ y5 (2) (r) { 9 ,u- -z,K) T= V, (r) . (6.1) 

(6.2) 



The noise correlator is found from (5.9) to be 



V dT V y a y a \ dT 



(4tt a ) z„ + — z 



rfr V 2v a v a s 2 . 
p. 



2T 2//3 
(4^V) 




dT' 



V 



VfiS 



(6.3) 



where we have substituted the vacuum Hadamard func- 
tion for the field-noise correlator (xx'} ■ The antisymmet- 
ric combination futvi vanishes, and hence, 



Mt)^(t')) =0, for 



t= 0. 



(6.4) 



The origin of this at first surprising result is not hard 
to find. First, the "z x" dependant terms vanish since 
the average acceleration is zero, but why should the other 
terms in the correlator vanish? From ( |6.3| ), we see that the 

second term in the particle-noise correlator involves z 
9\]G H (z, z'). But the gradient of the vacuum Hadamard 

kernel (which is a function of a) satisfies d^G H (a) oc y^, 
and is therefore always in the direction of the spacctime 
vector connecting z (r) and z (r') . However, for the iner- 
tial particle, y^ oc v^, and thus the antisymmetrized term, 
z\u d u ]G H (a) oc V[nV u ]G H , always vanishes. It therefore 
follows that a scalar field with the coupling that we have 
assumed does not induce stochastic fluctuations in a free 
particles trajectory. This result is actually somewhat a 
special case for scalar fields, because in Paper III || we 
show that the electromagnetic field does induce stochastic 
fluctuations in a free particle. In the case of the electro- 
magnetic field there is more freedom in FfjjJ* — d^A^ 
(e.g., is a vector field; d^ip is the gradient of a scalar 
field) and so one does not find the cancellation that oc- 
curs above. For the case of a scalar field, there are noise- 
fluctuations when the field is not in the vacuum state (e.g. 
a thermal quantum field), and/or when the particle is sub- 
ject to external forces that make its average acceleration 
non-zero (e.g. an accelerated particle). 

We therefore find that the free particle fluctuations (in 
this special case) do not obey a Langevin equation. The 
equations of motion for z^ are 



d 2 z ll (r) e 



dr 2 



(6.5) 



Eq. ( |6.lD has the unique solution 
'K (r) = 



(6.6) 



after recalling that g (2) (0) = ^'f (0) = 0. That is, the 
same initial time behavior that gives unique (runaway- 
free) solutions for z (the semiclassical solution) apply to z 
(the stochastic fluctuations). 

Let us comment on the difference of our result from 
that of [Q. First, the dissipation term in our equations 
of motion is relativistically invariant for motion in the vac- 
uum, and vanishes for inertial motion. The equations of 
motion in |46| ] are not; the authors concluded that a par- 
ticle moving through the scalar vacuum will experience a 
dissipation force proportional to its velocity in direct con- 
tradiction with experience. This result comes from an in- 
correct treatment of the retarded Green's function in 1+1 
spacetime dimensions- which is quite different in behavior 
than the Green's function in 3+1 dimensions. We avoid 
this problem by working in the more physical three spatial 
dimensions. We furthermore find that a free particle in the 
scalar vacuum does not experience noise (at least at the 
one- loop order of our derivation- see Paper I) . This result 
is somewhat special owing to the particular symmetries of 
the scalar theory. In Paper III, we find that charged parti- 
cles experience noise in the electromagnetic field vacuum. 
We have not addressed in depth the degree of decoherence 
of the particle motion, and how well the Langevin equa- 
tion approximates the quantum expectation values for the 
particle motion is still an open question that will depend 
on the particular parameters (e.g., mass, charge, state of 
the field) of the model. 

Our result also differs from that of Ford and O'Connell, 
who propose that the field-cutoff take a special value with 
the consequence that a free electron experiences fluctua- 
tions without dissipation p8fl . We have already argued 
why this assumption is unnecessary as a consequence of 
recognizing QED as an effective theory. We have also em- 
phasized the distinction between radiation reaction and 
dissipation. Hence, we find that a free (or for that mat- 
ter a uniformly accelerated) particle has vanishing aver- 
age radiation reaction, but when there exists a fluctuating 
force inducing deviations from the free (or uniformly ac- 
celerated) trajectory, there is then a balancing dissipation 
backreaction force. What is a little unusual about the 
scalar-field-free-particle case is that the "overly simple" 
monopolc coupling between the scalar field and particle 
does not lead to quantum fluctuation induced noise when 
the expectation value of the particles acceleration vanishes 
(i.e. for a free particle). In paper III [Q], treating the elec- 
tromagnetic (vector) field, there is quantum field induced 
noise and dissipation in the particle trajectory, for both 
free and accelerated motion. 
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VII. DISCUSSION 

We now summarize the main results of this paper, what 
follows from here, and point out areas for potential appli- 
cations of theoretical and practical values. 

Perhaps the most significant difference between this and 
earlier work in terms of approach is the use of a particle- 
centric and initial value (causal) formulation of relativistic 
quantum field theory, in terms of world-line quantization 
and influence functional formalisms, with focus on the 
coarse-grained and stochastic effective actions and their 
derived stochastic equations of motion. This is a gen- 
eral approach whose range of applicability extends from 
the full quantum to the classical regime, and should not 
be viewed as an approximation scheme valid only for the 
semiclassical. There exists a great variety of physical prob- 
lems where a particle-centric formulation is more adept 
than a field-centric formulation. The initial value formu- 
lation with full backreaction ensures the self-consistency 
and causal behavior of the equations of motion for the 
semiclassical limit of particles in QED whose solutions are 
pathology free- the ALD-Langevin equations. 

For further development, in Paper III we shall extend 
these results to spinless particles moving in the quantum 
electromagnetic field, where we need to deal with the is- 
sues of gauge invariance. We find that the electromagnetic 
field is a richer source of noise than the simpler scalar field 
treated here. The ALD-Langevin equations for charged 
particle motion in the quantum electromagnetic field de- 
rived in Paper III is of particular importance for quan- 
tum beam dynamics and heavy-ion physics. Applying the 
ALD-Langevin equations to a charged particle moving in a 
constant background electric field (with the quantum field 
in the vacuum state) , we show that its motion experiences 
stochastic fluctuations that are identical to those experi- 
enced by a free particle (i.e., where there is no background 
electric or magnetic field) moving in a thermal quantum 
field. This is an example of the Unruh effect fl47|| . The 
same effect exists for the scalar field example treated in 
this work, but we leave the essentially identical derivation 
for the more physically relevant case of QED presented in 

In a second series of papers 134] , we use the same con- 
ceptual framework and methodology but go beyond the 
semiclassical and stochastic regimes to incorporate the full 
range of quantum phenomena, addressing questions such 
as the role of dissipation and correlation in charged parti- 
cle pair-creation, and other quantum relativistic processes. 

In conclusion, our approach is appropriate for any situa- 
tion where particle motion (as opposed to field properties, 
say) is the center of attention. Our methods can be applied 
to relativistic charged particle motion not only in charged 
particle beams as in accelerators, and free-electron-lasers 
as in ion-optics, but also in strong fields such as parti- 
cles moving in matter (crystals) or in plasma media as in 
astrophysical contexts p8[. 



Self-consistent treatment of backreaction from quantum 
field activities on a charged particle is an essential yet often 
neglected factor in many problems dealing with charged 
particle motion. First principles microscopic approach by 
means of covariant quantum field theory leading to micro- 
scopic stochastic equation is an effective remedy to this 
deficiency. 

Theoretically, by interpreting a particle's spacetime 
properties (e.g. its location in time and space) as ef- 
fectively microscopic 'clocks and position markers' (e.g. 
events) we see that a quantum-covariant description of 
particle 'worldline' motion is not just particle dynamics, 
it is the quantum dynamics of spacetime events. For ex- 
ample, the "time" of the particle is treated as quantum 
variable on equal footing with position, and hence ques- 
tion like localization apply not just to localization in space, 
but localization in time. It would be natural to apply 
this approach to subjects like 'time of flight', that have 
become popular, and questions regarding decoherence of 
time (variables), together with space variables (e.g. the 
emergence of (local) time for particles). 
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FIG. 1. The trajectory and the lightcone of a particle at 
2 (r) . The radiation reaction force on a particle at z (r) de- 
pends on the particle trajectory in the interior of it's past 
lightcone. For a massless field, the radiation reaction force 
is local except that the regulated Green's function Ga smears 
out the radiation reaction over the past trajectory expanded in 
a Taylor series around z (r),as in Eq. (3.8). 




1 2 3 

FIG. 3. The time-dependent coefficient gi (r) that deter- 
mines the radiation reaction (RR) plotted against the ellapsed 
proper-time, r = r — n, since the intial time ti at which the 
factorized intial state of the particle plus field is specified. The 
radiation reaction vanishes at r = 0, but quickly builds to 
the asymptotic value familiar from the Abraham-Lorentz-Dirac 
equation on a timescale I/A. 



Proper-time (units of 1 /A) 




FIG. 2. The time-dependent renormalized mass m (r) of the 
particle plotted against the proper-time, r = r — n, ellapsed 
since the intial time ti at which the factorized intial state of 
the particle plus field is specified. The actual mass-shift de- 
pends on the cutoff A; the vertical mass units are arbitrary 
since they depend on the unkown particle bare mass mo. The 
renormalization time-scale is r ren ~ 1/A. 
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